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I INTRODUCTION 


When analyzing or designing a control system, it 1s important to study the effects 
On overall system performance of varying one or more parameters (mass, inertia, gain, 
resistance, etc). It is equally important to determune whether a desired dynanuc behavior 
can be achieved with any set of values for the parameters ... if not, redesign 1s indicated. 

The analvsis and synthesis of linear feedback control systems, or the compensation 
of same, can be realized by parameter plane methods. The parameter plane method, 
which works well for two variable parameters and which may be extended to three or 
more paranieters, is purely algebraic, and the resulting plots are valuable aids to analysis. 

Furthermore, the performance index cost curves in terms of integral-square error 
and other cost functions can be superimposed on the parameter plane graphically. That 
is very useful for analysis and for designing a control system when its characteristic 
equation has two parameters. Since the coefficients of the characteristic polynomial are 
deternuned by the system parameters, it follows that some relationship exists between 
the value of anv parameter and the value of the characteristic roots. The integral-square 
error as a performance index is also determined by the system parameters. 

Let us examune the choice of the integral-square error as a performance index. In 
niany applications the performance of feedback control systems can be considered either 
satisfactory if the errors are below a certain linut, or unsatisfactory, if the errors are 
above a certain linut. The hvpothesis that the configuration is fixed refers to the form 
of the system; it is assumed that one or more parameters can be adjusted. Given the 
input and ideal output signals as transient functions of time and given the configuration 
of the control system, find the values that the free parameters should have in order to 
nunimize the integral-square error. 

One approach to solving this problem starts with the Fourier transformis of the input 
and ideal output signals. The integral-square error is then easily formulated by means 
of Parseval’s theorem. If the transform of the error is a rational function of the complex 
frequency, the integral-square error can be easily evaluated by means of a _ table of in- 
tegral, which is in Appendix A. When the integral-square error is evaluated, it appears 
as a function of the free parameters. 

The object of this thesis, then, is to develop a user friendly interactive computer 


program to plot families of performance index cost curves. A secondary objective is to 


demonstrate the use of these curves in the design of the system. The design procedure 
starts with the performance index cost curves superimposed on parameter plane curves. 
The result of the analvsis is a choice of values for the parameters such that the svstem 


performance is the best available compromise between transient performance and cost. 


~ 


Il. BACKGROUND THEORY AND EQUATIONS 


A. FUNDAMENTAL CONCEPT 

Classical control theory, historically speaking, has been applied to the design prob- 
lem with the assistance of graphical presentation and trial and error methods. Such 
methods have been quite successful in the development of good control systems, but do 
not answer the question “Is this the best system possible.” The engineer and mathe- 
matician have always looked for the “best” contro] system. Early methods used “stand- 
ard forms.” The concept of a “performance index” was invoked, also at an early point, 
for, if the chosen index could be satisfied, one could say that the resulting system was 
fe@ecatie best’ as measured by the index. 

Let us take the performance index that is the integral-square value of the error be- 


tween the ideal output and the actual output. A simple feedback control svstem with 


adjustable parameters ts assumed as presented in Figure 1. 





Tigure Il. Block diagram of feedback control system. 


The performance index in terms of integral-square error 1s defined as: 


PI =|F dt (2.1) 


Where Pl> periormancesindex 


E ; input - actual output 


and the criterion was to design the system such that this index was minimized. For lin- 
ear systems this problem can be solved by applving Parseval’s theorem [Ref. I: p. 43] to 
obtain explicit relationships for system parameters. One should note that a system de- 
signed to minimize this performance index has a step response with considerable over- 
shoot and oscillation, and this may not be acceptable for some applications, but where 
such dvnamic characteristics are suitable the design is straightforward. Other perform- 


ance jndices which have been found useful are 


PI = {Ic td (2.2) 


PI = {le a 2.3) 


Both of these are nonlinear and do not have analytic solutions. Designs using such 
criteria Were accomplished using simulation. The solution of an optimal control problem 
is the result of minimizing the chosen cost function and is the design which best satisfies 


the criteria defined by that cost function. 


B. PARSEVAL’S THEOREM 
Suppose we wish to know the area beneath the product of two functions x,(z) and 
x,(t) of the infinite interval of the time from minus infinity to plus infinity. Let this area 


be denoted by the symbol I (which stands for integral), thus 


+00 
[= | x(t) x4(0) de (2.4) 


—oO 


Suppose, in addition, that the time functions are Fourier transformable and have trans- 
hommisex(S) aid x,(5), respectively. We wish to express the integral I directly in terms 
of these transforms in order to obviate inverse transforming x,(s) and .x,(s) when these 


are given rather than the time functions themselves. To do this we observe that 





+joo is . 
@ x4(s) ds >) 


Substitution of this value of x,(z) into Equation (2.4) vields 


a 2 1 +joo ee 
1=| x, (4) +] e@ X3(s) ds | dt (2.6) 


—co J) 5-6 


Let us now interchange the order of the integrations so that we integrate with respect 


to time first 





—joo —so 


SV fica eon, 
[= | X>(s) as| e x,(t) at (257) 


250] 


By the direct Fourier transform we know that 


x,(s) = [re X,(2) dt (2.8) 


—20 


Thus the integral with respect to time on the right side of Equation (2.7) mav be evalu- 


ated as 
+00 
| on x(t) dt = x,( —s) (225) 


This permits us to write Equation (2.7) as 


l +jo0 
I= | x,( —S) xo(s) ds (2.10) 


27j 





JOO 


This is the desired expression for the integral I. 


An important particular case of the above result occur when the two time functions 


involved are identical so that the integral I becomes 
FOO; oF 
I -| Sena: (220) 
—oO 


We call this integral the integral-square value of the time function .x(1). 


By Equation (2.10) we know this can be expressed as 





l +/joo 

l= >— | x( —s) x(s) ds (22) 
Lie tae 

Ji2,2) 

This result 1s Known as Parseval’s theorem. This theorem represents a very convenient 


way of expressing the integral-square value of a time function in terms of its transform. 


C. MINIMIZATION PROCEDURE 
n the basis of Parseval’s theorein we are now in a position to formulate the pro- 
cedure for minimizing the integral-square error for transient input signals by adjustment 
Ol the Tree par anmevens: 
This procedure involves four steps. The first step is to express the Fourier transform 
of the error as a function of the complex frequency s. This function will involve the free 
parameters of the system as unknown coefficients. The second step is to express the 


integral-square error I, in terms of the error transform E(s) by Parsevals theorenmenes 


+joo 
-_ 7 | E( —s) E(s) ds (2.13) 


—joo 


At this stage, provided E(s) is a rational function, the integral-square error will appear 


in the form 





{ c( —s) c(s) 


~ 2nj d( —s) d(s) a 


e 
Where c(s) and a(s) are polynomials in s. The third step is to evaluate the integral. 
Fortunately, definite integrals of this form have been evaluated in terms of the coefh- 
clients appearing in the polvnonuals. A table of these integrals along with an explanation 


of how they were evaluated is given in Appendix A. Upon carrying out the integration 


memimave the intesral-square error expressed as a function of the free parameters 


p,through p, Symbolically the integral-square error at this stage 1s expressed as 
1(P1, Po, *°* + Px) 


The fourth step of the minimization procedure is to adjust the values of the free pa- 
rameters in such a way as to minimize the integral-square error. Formally this can be 
done in the standard manner of equating the partial derivatives of the integral-square 
error to zero and solving the resultant set of equations for the values of the parameters. 


In svmbols this means solving the set of k equations of the form 


P 
Gt, 





CPr 
Unfortunately in many practical problems this formal procedure for finding the param- 
eter values that minimize the integral-square error leads to sets of non-linear equations 
in the parameters which are so complex that straightforward methods of solution do not 
est. 

If the equations for the parameters obtained bv setting the partial derivatives of the 
integral-square error with respect to the parameters equal to zero are too complex to 
be solved by analytical methods, we may resort to trial-and-error or successive approxt- 
mation techniques. However, if such methods are necessary to solve the equations for 
the parameters, it is probably best to go back to the expression for the integral-square 
error itself and plot this as a function of one of the parameters, holding the other pa- 
rameters fixed. By using a sufficient number of fixed values for these other parameters 
and making a sufficient number of plots, values can be found for all the parameters that 
define the minimum value of the integral-square error closely enough for engineering 
purposes. By working directly with the integral-square error expression, We avoid taking 
the partial derivative of this expression with respect to each of the free parameters. 


Taking these derivatives mav be a very difficult and time-consuming chore. 


Hl. PERFORMANCE INDEX CURVE EVALUATION 


A. BACKGROUND 
l. Overview 

The performance index program Is an interactive personal computer program 
whenever possible. Before using this program, we have to generate a datafile containing 
a. fPand PI. This file will enable us to plot families of performance index cost curves. 
There currently exists, at Naval Postgraduate School, a DSL(Dynamic Simulation Lan- 
guage) [Ref. 2] program residing on the school’s IBM 370 mainframe. By using the DSL 
program, we can easily obtain the above mentioned data file, and then, down load this 
ile from the mainframe computer to floppy diskette. This data, now on floppy diskette, 
is to be used, along with the computer program developed in this research paper, to plot 
performance index curves. With this program, the user can examine particular values of 
the performance index curves. Options, such as printer type and size of output, can be 
selected Os the wer 

2. Software / Hardware 

The program language used in coding this program is Microsoft FORTRAN77 
V3.20. The Plotworks PLOTS88 graphics librarv [Ref. 3] is used to generate output plots. 
The plotting of performance index curves conducted for this thesis were performed on 
an IBM-AT or IBM compatible $0286 machine. Due to the universality of FORTRAN 
coding, this program could be implemented on any machine capable of being pro- 
grammed in FORTRAN. The source code of generating curves of constant perform- 
ance index is listed in Appendix B. 

Before running the program, ANSI.SYS should be incorporated in the personal 
computers CONFIG.SYS file prior to running the performance index program. 
Through the ANSI.SYS device driver, system calls to clear the screen and position the 
cursor are enabled. If not, screen readability suffers although the program ts still fully 


functional. More will be said on this in the ANSI module section of this chapter. 


>, 


Se 


Figure 


Iwo Dimensional Contouring 


The control grid is constructed as follows Figure 2. 


tshech yheehs 


g 
& 
E 
= 
z 


tdi 


MVvsizZt 
pony 


2. he contour grid two-dimensional real array z from Ref.3. 





nusize 


where +» = dimensional size of real array z. 


(nX.nv) = maximum index of area used. 

(xlow.ylow) = lower left corner of the grid area in the world units. 
(xhigh.vhigh) = upper nght corner of the grid area in the world units. 
(nxsize.nysize) = grid interval in world units. 

dx = (xhigh-xtow),/(nx-1) 

dv = (vhigh-ylow)’(nv-1) 


a 


Note that only the control points in the range (i= 1,.....,.nx ;]= l.,......ny) are used 


in the array z. By convention, the first index is referred to as the row index and the sec- 


ond index is referred to as the column index. 


The following sequence produces a two-dimensional contour map: 


A genid constructed over the sunace Hicuress 

At each control point of the grid, the z ordinate is determined. Using Laplace or 
Spline interpolation or a combination of the two. curves of constant z are con- 
structed bv PLOTS8S graphics hbrarv. To do that, we can use the subroutine named 
ZGRID which 1s the PLOTSS hbrary, this subroutine interpolates a set of arbi- 
trarily spaced data points xp(X), yp(k), k =1, mpts onto a rectangular grid z(i, /) of 
arbitrary S1Ze. 


The grid control points may be smoothed. The ZSMTH subroutine which is the 
PLOTSS library smooths the gridded control points by appiving Laplacian 
smoothing. 


Potions of the contour may be removed to make room for annotation. 


The contours are drawn. Each set of 4 grid points can be subdivided further into 

4+ or 16 subgrids for more precise contouring. In this program, the set of 4 grid 
points 1s subdivided into 16 subgrids. Variable contour line tvpes and widths can 
be drawn. To make drawing more precisely, we can use the subroutine named 
ZCSEG which is the PLOTSS library. This subroutine converts the grid of contro} 
points into a set of contours and draws the contours. In addition, the capabihtv 
of dividing each grid square into 1, 4, 16 subsections using cubic polvnonual in- 
terpolation. This step may be needed if contours are close together or if the line 
quality of contour need improvement. 


The size of the contour grid is defined in the main program. For most applica- 


tions, a size of 50 x 50 control points is sufficient. The maximum number of control 


points is limited by the available space in the data segment of the program. A maximum 


of the 51 control levels can be drawn. 


10 


EOE _-t 


B. THE PROGRAM 
The performance index program given in Appendix B ttself consists of main program 
and six subroutines which make up the performance index package. Figure 3 1s a flow 


chart illustrating the organization. 


START 


a 


processing 


display 
monitor 


Hardcopy ? 
No 





Figure 3. Performance Index Program Flow Chart. 


When the program ts initiated, the following questions are presented: 
1. Ilow many constant PI curves # = = = > 


e The user has to enter an integer value which is must be less than 51. 


| 


2. Enter PI valucs, lwet, ldig # = = = > 


e Table 1 give us proper values. 


Table 1. PARAMETER TABLE 


parameter value Description 


lwet Integer [tomer the lmnve tvpe for cach PI! curve the Pleune 
is drawn with a solid line. See PLOTS8 
Vidi cle 53 | 


Idit omit decimal point. 
mteger = ont contour annotation for this level. 





Ee tOese the PI curves ts drawn with a dashed line. 


Idit integer > 0 the annotation control information for each 
Pl. number of digits after decimal point. 
< -2 


3. Enter number of data file points = = = > 


e The user has to enter total number of data points. 


I. Subroutine SIZE 
The subroutine SIZE Figure 4 is displayed. This subroutine permits adjustinent 
of plot size to both the screen and the printer. This option was utilzed to scale plot 


outputs for melusion in this thesis. 


Output plots are currently sized to fill the console screen. 


Would you like to adjust the plot size? 


LA + 9 9 


Enter “y or n =>) y 


A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 
It will be doubled with a factor entry of 1.8. 


Enter’factor size ==> 0.7 





Figure 4. The plot size subroutine. 
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2. Subroutine LDFILE 


iersuiproutme LIDTiLE (Figure 5) 1s displaved on the screen. 


File name should not exceed 8 characters 
File extension should not exceed 3 characters 


What is the file name (fn) and extension (ext)? 


Enter in form fn.ext (e.g. PARAPLAN. INP) 





Figure 5. =The data load subroutine. 


Subroutine LDFILE permits program input of a pre-existing data file. The user 
is queried as to file name (not to exceed eight characters) and file extension (three char- 
acters or less). The code checks for existence of the file within the working subdirectory. 
[fit does not exist, a message 1s returned stating such and offering the option of entering 
another file name again. Entry of file extension is optional. In fact, a file can be labeled 
maithens little as one character or number. 

We note that the data should be formatted tn order to read the data file suc- 
Secsimily by: Using the {DFILE. A data file could be created aligning in the format 
identified by the subroutine LDFILE source code listed in Appendix B. 

3. Subroutine NMONPRT 

The largest of the subroutines 1s subroutine MONPRT. A number of smaller 
subroutines are associated with it. Subroutine MONPRT allows selection of the graphics 
Output device, be it the monitor or a wide variety of printers. Graphics output defaults 
to the monitor when the program is first executed. Graphics output can always be 
dunmiped from the screen to a printer, but there may be times when direct output to a 
printer is desired. All NPS printers are included in the printer output options. In addi- 
tion, One Option permits the user to directly enter the values for [OPORT and MODEL, 
as outlined in the PLOT88 manual. [Ref. 3: p. 17] 


The first menu displayed when subroutine MONPRT is called is the 


printer/output menu Figure 6. 


Epson FX-80, All 

Epson FX-100, All 

Epson MX-100, All 

Epson RX-80, All . 

Epson MX-80 & IBM Printer 

HP 7470A Graphics Plotter 

HP 7475A Graphics Plotter 

HP 758xB Series Plotters 

HP 2686A Laser Jet 

HEFL HEHEHE HEE EHEEE HEHEHE EE EEE EEE EE EEE EEE FFE THEE EE EEE eee ttt eeetet 
Graphics Monitor (default) 

HARDWARE Interface Menu 

Input PLOT88 Values for IOPORT and MODEL 
RETURN to LOAD Data 


i 
2 
3 
4 
5 
6 
7 
8 
9 


Enter integer number for selection == 





Figure 6. The printer/output menu. 


Selection of a printer in this menu automatically directs graphics output to the most 
commonly associated output port, either parallel or serial, for that particular device. 
Hlowever, usually there are multiple parallel or serial ports to which a printer can be at- 
tached. The default printer output setting mav not direct plots to the appropriate port. 
To provide maximum output flexibility, the printer/output menu allows selection of a 


specific output port through access to another menu. 


l4 


Printer Port 

Printer Port 

Printer Port 

Serial Port 

Serial Port 

FEEEEEEE EEF EE EEE EEE EEE EEE EEE EEE EEE EE EE EERE EEE EEE EEE EEE EEE EH HTH 
RETURN to LOAD Data 





Figure 7. Tlie Hardware Interface Menu. 


The output port selection option is called the HARDWARE INTERFACE 
MENU Figure 7. This menu is part of a subordinate MONPRT subroutine titled 
PORT. The option allow selection of one of three parallel ports (LP T1 - LPT3) or either 
of two serial ports (COMI - CONM2). 

Ifa sertal port 1s manually selected in the HARDWARE INTERFACE MENU, 
an associated data transfer (baud) rate must also be assigned. Immediately after se- 
lection of serial port options 4 or 5, the BAUD (data transfer) RATE MENU 1s dis- 
plaved Figure 8. Transfer of graphics data over a single line is usuallv time consuming, 
therefore the highest printer capable transfer rate available (9600 baud) is normally the 
best selection. A slower transfer rate would only be indicated if improper graphics out- 


put 1s generated. 


Enter integer number for selection == 





Figure 8. Baud (data transfer) Rate Nenu. 


Bit by bit transfer of data can also have a check sum (parity) associated with it 
to provide an interval check of correct data transfer. Following selection of a data 
transfer rate, the PARITY MENU Figure 9 is displaved. The user has the optionser 


choosing odd,even or no parity. 


NO Parity 
EVEN Parity 
ODD Parity 


Enter integer number for selection == 





Figure 9. ~—- Parity Menu. 


The PLOTS8 Software Library Reference Manual [Ref. 3: p. 18] contains a 


number of tables which provide common output settings for many more printers than 
canbe selected inthe PRIN BER: OU Wie esta iu 


C. ANSI MODULE 

The ANSI module contains two subroutines. Subroutine CLR clears the monitor 
screen. It is used primarily as a precursor to displaying a menu or positioning queries 
or informational text. Subroutine CUP positions the cursor immediately after a request 
for response. This subroutine serves not only to give a pleasing screen display but also 
provide the user with an inumediate reference as to where he is the program and if input 
is required of the user. Both routines can be found in the PLOTS88 Reference Manual 
[Ref. 3]. The two subroutines which comprise the ANSI module make use of ANSI es- 
cape sequences. These sequences are a series of characters that can be used to define 
functions to MS-DOS. In addition to clearing the screen and positioning, moving, and 
saving the cursor’s position, ANSI escape sequences can control screen graphics and 


reassign kev definitions. 


In order to utilize the escape sequences, the microprocessor must first be able to 
fieroret the commands. It does this through a systems file entitled ANSI.SYS. 
Pewoleo oS must be defined as a device in the CONFIG.SYS file to enable proper inter- 
pretation of ANSI escape sequences (e.g., device= \ for \ ansi.svs where \ for \ is the 
Mmmpetonasi.svs). If ANSI.SYS is not resident in the CONFIG.SYS file, the perform- 
ance index program will not clear the screen or position the cursor when so directed. 
Instead, the indicated escape code sequence will be printed on the monitor whenever a 
call is made to CLR or CUP. Although program readability 1s degraded, the perform- 


ance index program will still function correctly. 


I 


IV. PARAMETER PLANE DESIGN -GRAPHICAL METHOD 


A. OVERVIEW 

The algebraic solutions have the disadvantage that a fixed value of zeta and omega 
must first be chosen to compute the alpha and beta terms. In some instances it is pos- 
sible to modify the remainder polynomial so as to ensure that the specified roots are 
dominant. However, it 1s not always possible to guarantee that roots placed at a speci- 
fied location can be made dominant. Thus, an exhaustive trial-and-error procedure may 
be required to achieve the best values for the various parameters. To avoid this trial- 
and-error analytical approach one can employ the graphica! solution. 

We have two programs available. One permits us to get constant zeta curves and 
omega curves.The second program gives us a families of constant performance index 
curves in the parameter plane. We can incorporate the results of two programs by 
superimposing the plotted performance index cost curves over the parameter plane 
curves on the same sheet. 

To do that. we perform a preliminary procedure. This procedure involves four steps. 
The first step is to plot the parameter plane curves for the specific problem. This curve 
permits us to choose ranges of two parameters that are to be used in the analysis and 
design. The second step is to obtain a data file containing « , 8 and performance index 
values by using the D.S.L program. We then copy this data file from the mainframe to 
a floppv diskette. The third step is to plot the families of performance index curves in 
the same range of parameter values on the two parameter plane. We then have two 
kinds of parameter plane curves. The fourth step of the analysis for a given specific 
problem is to duplicate these curves on the separate transparent sheets and then put 
them together on the same sheet. Recapitulation of the graphical analysis path is shown 


in Figure 10. 
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Figure {0.9 Graphical analysis path. 





B. GRAPITCAL SOLUTIONS 

In this section, three different kinds of control problems will be demonstrated. One 
of these is a cascade compensator system. The other one is a velocity and acceleration 
fecdback compensator system. The above two problems will be evaluated with three 
mds of performance index constraints, ic, Pl=fEdr , Pl=fl£| td: and 
a f | E| dt. Vhe third problem is a steering control for optimized autopilots. The 


controller parameters are adjusted to muimuiniuze the cost function. 


ip 
P= + (i? + 5?) at 
0 


where 0 = rudder angle 
w= yaw error 


A = a weighting factor 
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1. Example 4-1 (cascade compensator) 
Figure 11 gives the block diagram of a high gain third order system which is to 


be damped with a cascade filter. [Ref 4: p. 319] 


Pl 














100 
5($+1)(s+10) 


Ee P (Ss + Z) 


Z (§ +p) 


Figure I]. The third order cascade compensated system. 


Eet 
eae 
yas 
B=P 
The characteristic equation for this system ts then: 
S* + (11 + B)S? + (10 + 11f)S? + (108 +100a)S + 1008 = 0 (4.1) 


Separating the coefficient terms. 


e Constant coellicient: S* + 1183+ 10S? 
© ALPIIA coefficient: 100S 


e BETA coefficient: S? + 11S? + 10S +100 


Enter each coeflicient array value and desired zeta and omega values into the computer 
program. igure 12 provides a look at the combined zeta and omega graph spanning a 


frequency range from 1.0 to 30, ¢ = 1.0, 0.5, 0.0 and w, = 5.0 
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Figure 12. 
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Combine zeta and omega curves for example 4.1 
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By inspection of the two paraimeter range, the alpha values vary from 0.0 to 140. 
and beta values also vary from 0.0 to 100. Now, within these ranges, simulation in 
D.S.L is used to evaluate the cost at a large number of points on the parameter plane. 
The D.S.L simulation program and its resulting data file containing « , f and three dif- 
ferent criterion cost values are given in Appendix C. 

Once we have obtained a data file, we download this fite from the mainframe 
computer to a floppy diskette. To run the performance index prograin, type in the exe- 
cutable file name, COSTFUN. The IBM microcomputers in the NPS Controls Lab are 
set up under the file management system IDIR. To execute the program under the this 
system, enter the appropriate subdirectory by using the up/down arrow keys to position 
the cursor adjacent to the subdirectory naine, then depress the ENTER key. Use this 
same procedure to select COSTFUN.EAE and push ENTER to run the prograniaine 
preparauion step for the plotting desired curves will appear. 


Enter proper values after each prompt as follows: 


Vv 


1.How many constant PI curves # 


2.Enter 
~Enter 
~-Enter 
-Enter 
~-Enter 
~-Enter 
-Enter 


jh 0M FO NM NO 


values 
values 
values 
values 
values 
values 
number of 


,»lwgt,ldig 
,»lwgt,ldig 
,»lwgt,ldig 
,»lwgt,ldig 
,lwgt,ldig 
,»lwgt,ldig 


data file points 


te te te te He Fe 


WvvvvvVv oO 
noooeo°oe°oe 


1 
2 
3 
4 
9) 
6 
t 


Next, adjust the plot size. 


—ean of oe oe oe eo 6 ee 2 oe ee ee ie 6 ie © 2 em oo ooo © es of eo eo ee © 8 ew ow © Oo 9 © of oe @ 6 oe © of ee ow ow oe & oe = 


Output plots are currently sized to fill the console screen. 


Would you like to adjust the plot size? 


Enter “y” or nn 2 == 5, 


A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 
It will be doubled with a factor entry Of wie. on 


Enter factor size ==> 0.7 


The last step is to enter a data file name. 


a eins 


File name should not exceed 8 characters 
File extension should not exceed 3 characters 


What is the file name (fn) and extension (ext)? 


Enter in form fn.ext (e.g. PARAPLAN.INP) ===> pil.dat 


After displaying the graph on the monitor, the program prompts you for a 
hardcopy of the graph. 


WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N) 


ee ee ee ee ee ee re ee a ee ee eee ee ed ee ee ee ee ee se ee ee Sh, SUL UL hl) 
FN I ee ee ee ee ee Oe ee Oe le le hm 


: 1 Epson FX-80, All : 
: Z Epson FX-100, All 
3 Epson MX-100, All 
4 : Epson RX-80, All 
5 Epson MX-80 & IBM Printer : 
: 6 HP 7470A Graphics Plotter 
7 : HP 7475A Graphics Plotter 
; 8 : HP 758xB Series Plotters : 
: 9 1 HP 2686A Laser Jet : 
FEES EHEE HEE EEE EEE EP EH HEHEHE HEHE HHH Heese tee tetettttetttt+tttetettt+ 
10 ; Graphics Monitor (default) : 
; 11 ' HARDWARE Interface Menu 
: 12 ; Input PLOT88 Values for IOPORT and MODEL 
' 99 RETURN to LOAD Data : 


bard ren ee ce ee ee ee Oe Te ee ee ee ee he he he he 


Enter integer number for selection ===> 5 


Reenter the data file name for the hard copy. Figure 13 is a performance index 
curve for constraint (ye dt. In like manner, the other two nonlinear criterion of per- 
formance index in terms of f|E| ¢dt and f|E]| dr can also be plotted by running the 
same program. Figures I4 and I5 are shown for desired constraint es tdt and 


JIE] dt, respectively. 
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Figure 13. The performance index curves for f £? dt. 
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Figure 14. 
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The performance index curves for {| £| ¢ dt. 
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Figure 15. 


00°00) Q0°08 00°09 00 'OR O0 0c 
B}ag 


The performance index curves for [| E| dt. 
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PemeUseconsider the two kinds of parameter plane curves, Figure 12 and 
Figure 13. By superimposing one over the other, the curves can produce the very in- 
teresting graph of Figure 16. In like manner, Figure 12 curves superimposed on the 
other cost contour are shown on Figure 17 and Figure 18. These graphs contain three 
kinds of constant families curves: € , w, and cost curves. At this point, from Figure 16, 


we could select arbitrary operating points. The selected ALPHA BETA pairs are: 
A point: « = 25, f = $0 
Bape o— 19505 — 60 
C point: o = 9.0, 6 = 23 
D point: « = 100, £ = 70 


Now, go back to the characteristic equation Eq. (4.1) and substitute in this equation 
each ALPHA;BETA pair and then, simulate the time response for each one. The tab- 
ulated simulation results are shown in Figure 19 and Table 2 on page 32. From the 
simulation results, we find that operating point D (« = 25, # = SO) had the fastest setting 
time, the lowest overshoot and the lowest cost index. The cost associated with this point 


are: 


[ea =0.11 
fiz Vale 
fre dt = 0.23 


And step response is given in Figure 19. In this case, a cascade filter was designed with 


the zero equal to 3.2 (= 80/25) and the pole equal to 80, using «= P’Z and f = P. 
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Figure 16. 
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Combine zeta, omega and cost curves for f E? dt. 
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Figure 17. 
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Combine zeta, omega and cost curves for [{ [E| t dt. 
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Figure 18. 
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Combine zeta, omega and cost curves for [| E| dt. 
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Figure 19. 
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Table 2. SINIULATION RESULITS. 
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The exact values, however, will be chosen to achieve a compromise between 
transient performance and cost, based on given system specifications. It is of interest 
to compare the performance index between the constraint f £? dt and the constraints 


f\E| edt and f|E| drat the same operating points. 


2. Example 4-2 (feedback compensator) 
This example deals with a plant incorporating both velocity (tachometer) and 
acceleration feedback. The block diagram of such a system is depicted in Figure 20. [Ref. 
d: p. 331] 





Figure 20. Velocity and acceleration feedback compensated system. 
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The characteristic equation for this system 1s: 


eres S0K2)S* + (24 SOKI)S £ 50 =—0 (4.2) 
eet 
Cool 
B= kK2 


Using the same procedure as the previous example, Figure 21 provides the combined 


zeta and omega graph spanning a frequency range from 1.0 to 30. 
€¢=0.45 
w, = 2.24 


This graph shows the two parameters range, the alpha values vary from 0.0 to 1.6 and 
beta values vary from 0.0 to 1.0. 

Now, we evaluate the cost at a large number of points on this range by using 
D.S.L. The D.S.L simulation program and its resultant data file containing « , 8 and 
three different cost values are given in Appendix C. Once we have the data file, we can 
produce the family of performance index curves. These are shown in Figure 22, 


meeure 23 and Figure 24, respectively. 
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Figure 21. 
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Constant zeta and omega curves for example 4.2 
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The performance index curves for f E* dt. 
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Figure 24. The performance index curves for {| E| dt. 
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Let us consider the parameter plane curves, Figure 25, Figure 26 and Figure 27. These 
were produced by superimposing Figure 21 Over each cost contour curve, Figures 
Figure 23 and Figure 24. 

Figure 25 consists of three families of curves: zeta, omega and cost curves. We 


could select arbitrary operating points given by: 
A point: « = 0.45, f = 0.18 
B point: « = 0.60, 6 = 0.40 
C point: « = 0.6, 6 = 0.067 
Dencink: 4 — Osa — 07 
E point: « = 0.4, 6 =0.125 


The operating point C (¢ = 0.6, 6 =0.067) permits us the lowest cost index. [he cos 


associated with this points are: 


fe dt =0.47 
fie tat = 0.25 


fre dt = 0.68 


However, settling time was too slow, ¢, = 0.6 sec. In this system, operating point C 
(a = 0.6, 8 = 0.067) determines the parameters Kl and K2, 
K1 = 0.6 


K2 = 0.067 
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Combine zeta, omega and cost curves for INEZ dt. 
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Figure 27. Combine zeta, omega and cost curves for {| E| dt. 
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Step response on exaniple 4.2 


Figure 28. 


Table 3. SIMULATION RESULTS. 


op- PI values 
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3.067 





Table 3 and Figure 28 show stimulation results for operating points A through E. The 
designer can then choose the operating point (a, #) depending on his desired specifica- 


trons. 
3. Example 4-3 (optimization of controller) 


As a final engineering example, consider the steering control for a surface ship. 


The block diagram of such a system shown in Figure 29. 


CONTROLLER 





Figure 29. Optimization of controller. 
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Given that the cost function is defined: 


J=+] (iw2 + 6 *) dt 


where 6 = rudder angle 
w= yaw error 
A = a weighting factor 
In this study, we used values given by R.E Reid [Ref. 5] for the SL-7 container 


ship. Reid uses a second order Nomoto model for the SL-7, 


_ i 
SG) ea 
and also uses a controller 
Ki(sT1l 4+ 1) 
OS) ere 


rp 


His results are given in Table 4. 


Table 4. REID RESULTS. 


- ae 


nos | K | Tin Sd 







Let us assume T1 and T2 are adjustable free parameters. 


6.G= 0.1084 0.4556(sT1 + 1) 
"we" §(90.36s+ 1) (sT2+1) 

l 

T2 

is characteristic equation for this svstem 1s: 


Let aa, —— 


90.36S° + (1 + 90.368)S* + (B +0.049a)S + 0.0498 = 0 (4.3) 
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The parameter plane curves examined within the frequency range from 0.01 to 0.1 
€ = 0.5, 1.0 window size; O< a < 35.,0< 6 < 1 Using the parameter plane program, we 
get the praph of Figure 30. 
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Figure 30. Constant zeta curves for example 4.3 
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In order to evaluate a family of cost curves, We have to obtain a data file on the same 
range of parameters, as figure 30 by simulation. The D.S.L simulation program can 


be implemented with ease by rearranging the block diagram in Figure 31. 


T1_(s+1/T1) — 
T2  (s+1/T2) elles OS 





Figure 31. Rearrangement of block diagram. 


The controller block becomes 


Ges 


(s + f) 





The D.S.L simulation program and data file for this example are given in APPENDIX 
C. Applying the previously discussed procedures for producing constant cost curves 


provides the graph of Figure 32. 
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Figure 32. The cost curves for - f TW? + 6) dt. 
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From Figure 32. we could select the operating points given by: 
A point:¢ =9, f=90.1 
B point: « = 20, § =9.2 
C point: « = 20, f = 0.4 


The operating point A («=9.0, 8=0.1) provides the minimum cost. 
Recall that 


R 
| 


S[- ap 


The two parameters values become 


T1 = 90.0, [2 = 10.0 


Comparison of Reid’s result and the simulation result show that they are very closelv 


matched. In addition, we simulated vaw response and rudder activity at the operating 


points A , B and C in Figure 33. Figure 34 compares the vaw response to course 


change for the ship at each operating point. Figure 35 compares the rudder activity. 


It is noted that operating point C is a faster controller but has a higher cost value. This 


cost is a higher rate of fuel consumption used to overcome the higher hydrodynamic 


drag forces. 
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Figure 33. 
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Figure 34. 


The yaw response ys time. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The graphical solution of the parameter plane technique provides the designer with 
a visual means of deducing how the dominant roots of the characteristic equation move 
about in the s-plane as two user-defined parameters defining compensator attributes are 
varied. This was illustrated in the example problems in Chapter 4. The graphical 
method is able to analyze easily, not only the linear constraint [£? di but also the non- 
linear constraint f|E| ¢de, f|E| de and any kind of cost function using the perform- 
ance index program. 

Constant cost value contours can be plotted in much the same presentation as the 
parameter plane plot. This procedure is a useful design tool which permits simultaneous 
analysis of both dynamic response and cost. 

As has been demonstrated, rapid solution of a wide range of linear automatic con- 
trol svstems problems can be achieved through use of the performance index program 


developed as part of this thesis. 


B. RECOMMENDATIONS 
Further work should be done to combine individual programs so that both param- 


eter plane and cost curves are produced by one consolidated program. 





APPENDIX A. TABULATED VALUES OF THE INTEGRAL FORNMI 


Table A.1 gives the value of 7, for values of » from | to 9 where 


_ 1 ft els)e(-s) 
n 2nj > aS d(s)d( —s) a) 
and 
c(s) = a +eestCy (A.1 — 2) 
Us) = d,S" +eeetdy (4.1 — 3) 
table A. 
¥ 
0 * 
a Baad,’ 
] Chdo ar cody 
“7% 2dod,d, 
F c5dydj + (c; — 2¢9C,)dya; + C6 dydy 
3 SE eS 


2dpa;( —day == d,d>) ; 


-_ ¢3( dp dy + dodydy) + (cf = 2eyes)dydyds + (cf — 2cocy)dydydy + col —dydy + dydyds) 
: Qed —dyd? — dd, + dyaydy) | 


‘Ss - [cam + (cf — 2eyeq)rm, + (cf — 2cyC3 + Zegeg)nty + (cy — Wocy)iy + coms J: 


a) = —dyd; == d,dy, 
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My — —dpd; + a,dz, 
1713 = as (aM, = d,my), 
do 
Ng = = (a, = dmb), and 


A; or Aya, =< d,m3 == dsm). 


l= = [ cgrtg + (ci — 2C3C5)m, + (c3 — 20¢5C4 + 2¢,C5)rty + (c3 — 2¢,¢3 + ogc, eam 
6 


Ta, Lei - 2€9Cy)M4Co sms |; 
where 


My = (dy, — dyn, + dys), 


allt, 
de 
Wi dyads cade idan ate 


713 = dyds ae d,d°d- ao d,d,ds, 


mn, = (dom, — dy + dgm)), 


ue 
dy 
Ms = Pe (ad), — dy, + dont), and 


As — d)(dms — dy, + dm). 


ia a 


a re nG = 20 NC 2 Chey ine Gh — 2c Cq)ms + come J; 
i. 


54 


[ cg my sh (cs — 2¢4C¢)m, + (ci — 2¢3C€5 + 2¢5C6)m, + (ee 205 Cpe Ofte 


2cyce)ms | 


where 
Be ig A A 
Wy = @ (dsm, — dam, + ams), 


= — (dd, — dods)” + (dod; — dyd,)(dod, — dy dg + dds — a3d,), 


m 


— 


My 


(dyad, — dd.) —dods + d,dy) + (dd; — a, d,)(ad,d, — a,dg), 


m3; = — (dod, — dds)” + (d)a, — d,dy)(dsd, — ds), 


My 3 (d,m3; — dymy + der), 
Cs = o (dm, cy dyn3 =F dems), 


Me = 7 (d5115 — dyn, + d,113), and 
0 
A, = ad; we a3INs + asm, — asim). 


| 2 2 2 z 
I, = BA, [ c5 mip == (C6 — 2C5C7)m, Ar (cs a INERE? == 2€3C7)m =P (c4 — 2C3Cs ap 2C5C6 = ee Nice, || 


as 2 2 2 
ap 2A, [ {c; eee ee eet Cee itl (Cy — 2CjCe te 2CpCa iia 4 (Cp — 20,65 )*it, + exrrale 
where 


my = (dod; — dyds)( —dydyd, + dodyds + 2d} dg) + (dydy — d5)(dods + d,d3) 
+ d\d3d3(doa; — d\d,) — dy de(dods — d,d;) 
+ ( —dyd, + dds — dyds)(dyd + dj) — dydo(dy de + 3dpdyds) 
SANG (nue ana ite Dances. 


Mm = (dod — dydy)(dyd; — dds — dydsd, + dydsd;) 
+ (dydg — dydy)( —dydyds + dods — dydydy + dydy) 
— dydsd;(dyds — dydy) + dj dg(dod; — dyd,), 
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3 = =a, (dds = dyd,)° + ( —dsd + d,d5)(dyd, ds = dyd; + d,a,d, == d: dy) 
4: dnd; (—adpde daar 


ms = (—dsdy + dedy)(2dydy dy — dydyds + dydyds — d\d,) 
+ (— dydy + dyd;)(dydyd, — dydyd, + dds) — dod, 


Ms, —! ra (dm, a dmg oe det — dgm)), 


Me — = (d511¢ = dyin, al den; = dgmy,), 


m= oa (do, — dyms + dein — dgn3), and 
Ag = dy(di ity — d3Img + dsttts — day). 


l 2 2 2 2 


l 
as 2A, (Cy — 2¢3C5 + 2€5C€, — 2c,C7 + 2CoCg)my + (cj = Fes oe 2¢,¢5_2cpea item 
cas 2. 94 23) ico 2 ) Hea 1; 
20g (C3 £0) C3 2CyC4)M¢ Cy Colo JING Coltig |, 


where 


nS = (dj, — dy, + dni, — dm), 


2 2 
My = Ay(A, Aq — AAg + A3Ag + A,ag + 204A — As — Asa, — as) 


3 
+ a, — A — A347 + Ayas + 2ayaz) — ay, 
2 2 
My = Ay(A3Ag + AyAg — Asdg + Aga; — A7Ag) + Ay( — a, Aq + Ayag + az) — aay 
4 2 2 
M3 = A\(G34 19 + AyQyq + A7dq — ag) + Ay( — aAy9 + 474g) — aya7, 


3 
Mg = Ay (Asay + 2a7a;9 — gag) + a,(a7, = 4Q19) — @, 


Ms = >a (d,m, — dym, + dm, — dgm,), 
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nN, = om (dims — dying + dm; — dginy), 


Mm, = =" (dmg — dynts + deity — agin), 


ie - (dy, — dyntg + dems — dym,), and 


Ag = dy(d,mg — dy, + dsimg — dams + doa), 


where 

a, = dad, — dyad, 
A, = a\d,— ads 
a, = a,d, — ad; 
ag = ads — dyd; 
as = dd, — d,d; 
a, = dasa, — dd; 
Az = a, dy — dodg 
Ag = ad, — aydy 
Ag = asd, — Aydy 


Qo = asa a dgdy. 


aii 


APPENDIX B. PERFORMANCE INDEX PROGRAM. 


This appendix contains a listing file of the FORTRAN source code for the per- 
formance index program. This program is used to plot the performance index curves for 
desired values. 

The program consists of main program and subroutines. The subroutine SIZE per- 
mits adjustment of plot size to both the screen and the printer. The subroutine LDFILE 
permits program input of a pre-existing data file. Subroutine MONPRT allows selection 


of the graphics output device, be it the monitor or a wide variety of printers. 


we ou a % %& % ARRAY DEFINITION * * % *& % % 


CHARACTER*1 CHAR,CHG 

REAL Z(51,51) 

REAL XP(202),YP(202),ZP(200) 
REAL ZPIJ(200) 

INTEGER KNXT( 200) , ITEXT(1) 
INTEGER J,K,SEL 

REAL ZLEV(50) 

INTEGER LDIG(50) , LWGT(50) 
COMMON/PLOTTR/ IOPORT,MODEL 


ve ode ve ve ve ve VARIABLE INITIAL VALUE * * * * se or 


DATA XMAX/11. 0/,YMAX/8. 0/ 

DATA XLPLOT/1.0/,YLPLOT/0. 5/,XHPLOT/8. 0/,YHPLOT/5. 5/ 
DATA CAY/5. 0/,NRNG/3/ 

DATA HGT/. 12/,NDIV/4/,NARC/8/ 

DATA NSM/4/ 


we ov ve ve ¥ % ENTER INPUT PARAMETERS * * * % % %* 


10 CONTINUE 
CALL CLR 
CAT GUPGIEG® 
MODEL=99 
IOPORT=99 
Wey Raines Pe) 
WRITE(**, 305) 
READ(*,*) NLEV 
IF (NLEV. LT. 0. OR. NLEV. GT. 50) GO TO 13 
DO 22 K=1,NLEV 
WRITE(**, 304) K | 
READ(*,**) ZLEV(K),LWGT(K) , LDIG(K) 
22 CONTINUE 
WRITE(**, 300) 
READ(*,*) NPTS 
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Wiki (%* , 297) 


ve we ye ve ve *& BEGIN PROCESSING * * vw w% % % 


ZS 


200 


FACT=0. 9 
CALL SIZE (FACT) 


CONTINUE 
eee iD ILECNPIS,XP,YP,2ZP) 


INITIALIZE PLOT88& 
CALL PLOTS(0, IOPORT , MODEL) 


DEFINE OUTPUT WINDOW 
CALL WINDOW(0O. ,0. ,XMAX, YMAX) 


REORIGIN THE ENTIRE DRAWING 
wean tO) 1.0,0.5,-3) 


XPMAX=XP( 1) 
YPMAX=YP( 1) 
XPMIN=XP(1) 
YPMIN=YP(1) 
DO 3 K=2,NPTS 
IF (XP(K).GT. XPMAX) THEN 
XPMAX=XP(K) 
END IF 
IF (YP(K).GT. YPMAX) THEN 
YPMAX=YP(K) 
END IF 
IF (XP(K). LT. XPMIN) THEN 
XPMIN=XP(K) 
END IF 
IF (YP(K). LT. YPMIN) THEN 
YPMIN=YP(K) 
END IF 
CONTINUE 
DO 4 J=1,NPTS 
XP(J)=XP(J)*XHPLOT/XPMAX 
YP(J)=YP( J)**YHPLOT/YPMAX 
XHIGH=XHPLOT 
YHIGH=YHPLOT 
XLOW=XPMIN*XHPLOT /XPMAX 
YLOW=YPMIN**YHPLOT / YPMAX 
CONTINUE 
ZERO GRID ARRAY Z 
mon200 I=1,51 
DO 200 J=1,51 
AG 5 )=0. 0 
CONTINUE 


GRID SIZE PARAMETERS 


NX=8 
NY=6 
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Ce 


580 


590 


900 


DX=( XPMAX-XPMIN) /FLOAT(NX-1) 
DY=( YPMAX-YPMIN) /FLOAT(NY-1) 


CALL FACTOR (FACT) 

CALL STAXIS (2134, 20n 150 1m—)) 

CALL AXIS(1.0,0.5, Alpha ,-5,-7.,0ee Ne De 

CALL AXIS(1.0,0.5, Beta’ ,4)=5. 0,905 vein on) 

CALL ZGRID(Z,51,51,NX,NY,XLOW, YLOW, XHIGH, YHIGH, XP, YP,ZP, 
NPTS , CAY,NRNG, ZPIJ, KNXT) 


POST THE INePUR eG inis: 


CALL ZPOST(XLOW,YLOW,XHIGH, YHIGH 
XP, YP,ZP,NPTS,2, (1x, £5. 2) 


,XLPLOT , YLPLOT , XHPLOT, YHPLOT, 
see) 


IRTN=590 
---- PUT A BOX AROUND THE MAP -------------- 


CALL RLON(XLELOnaV Le lOme. 
CALL PLOT(XHPLOT, YLPLOT, 2) 
CALL PLOT(XHPLOT , YHPLOT, 2) 
CALL PLOT(XLPLOT, YHPLOT, 2) 
CALL PLOT(XLPLOT, YLPLOT, 2) 


TPCIRIN. £05200) "GOlOs200 


CALL PLOT(1.0,0.5,-999) 

CALL PEOUGIAO. O52) 

CALL FACTOR (FACT) 

CALL “STARS (41a 20 O ieen 

CALL AXIS(1.0,0.5, Alpha’ ,-5,-7.40mexPMino) 
CALL AXIS(1.0,0.5, Beta (54.55.0420. ,YeEMIN, Uae 


-- SMOOTH THE DATA BEFORE DRAWING 325 ooo oe 


CALL ZSMTH(Z,51,51,NX,NY,NSM) 

CALL FACTOR (FACT) 

GALL STAMUS (13.7) 2On ss Ore) 

CALL AXIS(1.0,0.5,'Alpha ,-5,-7. ,0. ,XPMIN, DX) 
CALL AXIS(1.0,0.5, Beta {4)35.0, 90m. wenn) 


CALL ZCSEG(Z,51,51,NX,NY,XLPLOT, YLPLOD AHPLOn yer 


1 ZLEV,LDIG, LWGT,NLEV,HGT,NDIV,NARC ) 


IRTN=900 
GO TO 580 
CATE PLOtt 0 One se) 


WRITE(*,*)' WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N)' 
READ(* ,401)CHAR 
IF( (CHAR. EQ. 'Y').OR. (CHAR. EQ. 'y')) THEN 
CALL MONPRT 
GOTO 25 
ENDIF 


WRITE(*,*)' WOULD YOU LIKE TO DO ANOTHER ? (Y or N)' 
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Creat, Ce) 


C) C2 Cag 0) C2 2 


READ(*,401) CHAR 
IF((CHAR. EQ. 'Y'). OR. (CHAR. EQ. 'y')) GOTO 10 


297 FORMAT( 1x) 


300 FORMAT( ' 3.Enter number of data file points ===>', ) 
304 FORMAT( ' 2.Enter PI values ,lwgt,ldig # ',i2,'===>', ) 
305 FORMAT( ' 1.How many constant PI curves # ===>', ) 


401 FORMAT(A) 


end of program 


STOP 
END 
SUBROUTINE CLR ve 
SUBROUTINE TO ERASE THE DISPLAY. ¥ 
CALL, CALL CLR ve 


SUBROUTINE CLR 

CHARACTER*1 C1,C2,C3,C4 

INTEGER*2 IC(4) 

MoU ALENCE (C1,1C(1)),(C2,1C(2)),(C3,IC(3)),(C4,1C(4)) 
DATA IC/16#1B, 16s5B, 16#32, 16#4A/ 


WRITE ESCAPE CODE TO DISPLAY 
WRITE(*,1) C1,C2,C3,C4 
1 FORMAT(1x,4al) 


RETURN 

END 

SUBROUTINE CUP * 
MOVE THE CURSOR TO LINE N, COLUMN M. : 
CALL, CALL CUP(N,M) si 


N - NEW CURSOR LINE ‘s 
M - NEW CURSOR COLUMN SS 


SUBROUTINE CUP (N,M) 

GreRACTER*1 €1,C2,C5,C8,LC(5) 

CHARACTER*5 CBUFF 

INTEGER**2 IC(4) 

MeO OnENGE (C1,1C(1)),(C2,1C(2)),(C5, 1C(3)) , (C8, 1C(4)), 
*( CBUFF , LC(1)) 

DATA IC/16##1B, 16#5B, 16%3B, 16#66/ 


L=10000+100*N+M 


WRITE(CBUFF ,2)L 
2  FORMAT(iS5) 
MEME (e 1)) G1, C2, LC(2),LC(3),C5,LC(4) ,LC(5),c8 
1 FORMAT(1x,8al, ) 
RETURN 
END 


=e FSF Fees esses Ss SSF SF US S| CUS CU Ss CU SET hCUh ET Uh TCU ET TC SBT CU ST hCU ET UhU EE UCU EE TCU Bw—_UCUMWwET_CU wT BT SFT CU ST SBT Se eT hUhr EE hUr EU Ur EE UhCUM ET TCU EO TCU ET CU ZT OCU SBT CU ET CU ET hCU ET TCU ET CC ESE TCU UT CUT hCUrPEwTCUrET TCU TCU ET TUCU ET ST CU ET CU ET CU ET TCU TCU ET hCU EU hCU wT ThCUM EU TCU VEV—TUhUM TU ETUCUh ET TCU OWT hCUhEhCU Sh 
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C SUBROUTINE LDFILE -- ALLOWS LOADING OF DATA FROM AN EXISTING FILE * 


c ae 
ie ee 

SUBROUTINE LDFILE (NPTS,XP,YP,ZP) 

REAL XP(202),YP(202) ,ZP( 200) 

INTEGER NPTS,I 

CHARACTER*12 FILEIN,ANS*1 

LOGICAL EXIST 

i JCARERGIER 

WRITE(*, 100) 

100 FORMAT( 1X, ////,10X, ' wn rr nn nr nn nn ne nn eee nen nee ee 
*------ ',/,10X, %** THIS IS THE ROUTINE TO LOAD DATA FROM A FILE ** 
weet 7 10x, ! ee ee oe So Oommen ‘ 
*//,10X,'FILE NAME SHOULD NOT EXCEED 8 CHARACTERS ' ,/, 

*10X,' FILE EXTENSION SHOULD NOT EXCEED 3 CHARACTERS’ ) 
WRITE(*,101) 
101 FORMAT(1X,////,10X,'WHAT IS THE FILE NAME (FN) AND EXTENSION (EXT) 
*? ' 7/,10X,'ENTER IN FORM FN. EXT (E.G. PARAPLAN. INP) ===> ') 
CAL CURCEon ow) 
READ(*,102) FILEIN 
102 FORMAT(A12) 
INQUIRE( FILE=FILEIN, EXIST=EXIST) 
IF (EXIST) THEN 
OPEN( 7, FILE=FILEIN, STATUS='OLD' , ACCESS=' SEQUENTIAL’ ) 
ELSE 
WRITE(*, 103) 
103 FORMAT( ///,10%. THAT FILE DOES NOT Pxis® sce 
* 'DO YOU WANT TO TRY AGAIN? (Y/N)') 
READ(*,'(A1)') ANS 
IFCCANS .EQ. "Y') SOR. (CANS = BO3es 9) a tHE. 
GOTO 1 
ELSE 
GOTO 2 
END IF 
END IF 
DO 190 I = 1,NPTS 
READC7 120) REC 1) gyre ec) 
120 -FORMATOIX, F10-2 3p. or Ome) 
190 CONTINUE 
CLOSE(7,STATUS=' KEEP’ ) 
2 RETURN 
END 
eee ee 
C cc 
C SUBROUTINE MONPRT -- ALLOWS DISPLAY ON MONITOR OR PRINTER ¥ 
ce me 
C SS oe tt OO Ooo mam ooccocoo coco oC OO ooo oS 


SUBROUTINE MONPRT 
INTEGER, 1OPORT SOD EL out 
COMMON/PLOTTR/ IOPORT,MODEL 


CALL CLR 


80 
81 
87 


82 
83 


84 


93 


WRITE(*,80) 
WRITE(* ,81) 

WRITE(*, 82) 

WRITE(*, 81) 

WRITE(*, 83) 

WRITE(*, 87) 

WRITE(*, 84) 

WRITE(*, 81) 

WRITE(*,93) 

FORMAT(1x,///,10x,' wn nn nn nnn ne =e 
Wee--------------- ,/,10x,' ',21x,'PRINTER/OUTPUT MENU' ,22x,' ') 


Co ) 

FORMAT( 10x, | t+++t+t+++4+t+4+++t ttt tHe ttt + ttt Ht tt ttt ttt ttt ttt 
w+t+t++t++4-+ | ) 

FORMAT(10x,' ',2x,'PRINTER NO.',2x,' ',20x,'PRINTER',19x,' ') 
FORMAT( 10x, ’ sp 1 ",2x,  ',4x,' Epson FX-80, All' 

same OX , 
oe LOX, 
* 25x, 
amy, 10X, 
cao 5xX, 
cay) LOX, 
* 26x, | 
oo > 10x, 
Plex , | 


3x, 2 Are a oc, i sjoiorn (Dome), 7 ILIL 

1 
t 
t 
t 
t 
t 

*,/,10x,' "43x," 6 ',2x,' ',4x,'HP 7470A Graphics Plotter’ 

1 
t 
t 
t 
1 
t 


oe 3 (52%;  ,4xqeepcon Pix Oman 1) 


eo aa A ",2x,' ',4x,' Epson RX-80, All’ 


1 


ie eam 5 ",2x,' ',4x,'Epson MX-80 & IBM Printer’ 


a7 x 
Boy LOX, 
comely x, | 
ey, LOX, 

* 18x, ' 
a, 4 LOX, 
*,24x,' ©) 

fieemat(lOx, | ,3x, 
ee 10 ',2x,' ',4x,'Graphics Monitor (default)' 


1 


ro Xa 7 2 Gk ey SheGraphiccamlotter 
hs oa 8 ",2x,  |',4x, HP 758xB Series Plotters’ 


ic ene 9 ",2x,' ',4x,'HP 2686A Laser Jet’ 


seh Mal ",2x,' ',4x, ‘HARDWARE Interface Menu’ 


mmelOx, ,3x,' 12 ',2x,' ',4x,'Input PLOT88 Values for IOPO 
*RT and MODEL’ 

ls | ! 

10x,’ ',3x,' 99 ee eee RE LURN Sto LOAD Mara, | 23x, 
FORMAT(1x,///,15x,'ENTER INTEGER NUMBER FOR SELECTION ===> ') 

CALL CUP( 24,56) 

READ(*,*) SEL 

Gsm. EQ. 1) THEN 


IOPORT = 0 
MODEL = 5 

Piste (Sh .EQ. 2) THEN 
IOPORT = 0 
MODEL = 15 

Diseelt (SEL .EQ. 3) THEN 
IOPORT = 0 
MODEL = 11 
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ELSE IF (SEL .EQ. 4) THEN 
IOPORT = 0 
MODEL = 1 

ELSE IF (SEL .EQ. 5) THEN 
IOPORT = 0 
MODEL = 1 

ELSE IF (SEL .EQ. 6) THEN 
IOPORT = 960 
MODEL = 20 

ELSE IF (SEL .EQ. 7) THEN 
IOPORT = 9600 
MODEL = 30 

ELSE IF (SEL .EQ. 8) THEN 
IOPORT = 9600 
MODEL = 80 

ELSE IF (SEL .EQ. 9) THEN 
IOPORT = 9600 
MODEL = 60 

ELSE IF (SEL .EQ. 11) THEN 
CALL PORT 

ELSE IF (SEL .EQ. 12) THEN 
CALMAGIER 
WRITE(*,100) 

100 FORMAT(1x,////,5x,' INPUT VALUE TO BE USED BY THE PLOT88 GRAPHICS 
* PACKAGE’ ,/,5x,'FOR IOPORT ==>',//) 
GCALEVCUPC] 42) 
READ(* ,**) IOPORT 
WRITE(*,110) 
110 FORMAT(1x,///,5x,'| INPUT VALUE TO BE USED BY THE PLOT88 GRAPHICS 
*PACKAGE' ,/,5x,'FOR MODEL ==>',//) 

CALE GUECI2 20) 
READ(* ,**) MODEL 


ELSE 
IOPORT = 99 
MODEL = 99 
ENDIF 


WRITE(*,98) IOPORT,MODEL 
98  FORMAT(//,10x,'PLOT88 WILL USE THESE VALUES TO OUTPUT GRAPHICS: ' 
*.//,10x, IOPORT = 14.10%) MODEL — sine) 
PAUSE 
RETURN 
END 


SUBROUTINE PORT 
INTEGER IOPORT, MODEL, SEL 
COMMON PLOTTR/ IOPORT,MODEL 


CALL CLR 


WRITE(**,79) 
79 FORMAT(1x,///,15x,'SELECTION OF A PRINTER IN THE PREVIOUS MENU AUT 
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80 


81 
87 


82 
83 


84 


23 


*OMATICALLY' ,/,10x,'SELECTS THE MOST COMMONLY ASSOCIATED OUTPUT POR 
*T, BE IT PARALLEL ',/,10x, 
*'(LPT) OR SERIAL, FOR THAT PARTICULAR DEVICE. IF THE PROGRAM GRAPH 
aes ,/, 10x, 
*' ARE NOT BEING OUTPUT CORRECTLY, USE THIS MENU TO PROPERLY ROUTE T 
*HE' ,/,10x, 
*'GRAPHICS DATA TO THE OUTPUT DEVICE. ',////) 

PAUSE 


GALL CLR 
WRITE(*, 80) 

WRITE(*,81) 

WRITE(* , 82) 

WRITE(* , 81) 

WRITE(* , 83) 

WRITE(* , 87) 

WRITE(* ,84) 

WRITE(*,81) 

WRITE(* ,93) 

FORMAT(1x,///,10X, occ ccc n ne nnn nnn nn nnn nnn nnn nnn nnn ee- 
i — — —— ",/,10x,' ',19x,'HARDWARE INTERFACE MENU' 
peeox, |) 

FORMAT( 10x, | ---- o-oo on nn nn nn nnn nn nnn nn nn nnn ne eo eee 
oa on oe s 

FORMAT( 10x, | HEEE EEE ETE T EHH EHH EH HH Ht Ett tt tt tt tt ttt tt tt tt 
KHt+tt+t+++++! ) 

FORMAT( 10x, ' 92x, SELECT NO. ' ,3x,'--1.,20x..  PORI 22x," ) 


FORMAT( 10x, ' ose 1 2x. Galo ber Erintersbort 
el oOX, 
*,/,10x, i oe 2 "2s. |; [oe Rt? Pinter Port: 
7 15X,.. 
eeeilox.' | ,3x,' 3 nox (oie heeeeinter Port 
5x, 
pmelOx,' ',3x,' 4, 2 eee Olicentale Pore 
Selox,: | 
ee ox,' |, 3x,’ 5 exes ae xem CON? (sexta pou 
* 16x, | ') 

FORMAT(10x,' ',3x, 
*' = 9g mou 14x, RETURN to LOAD Data.,13x, ») 


FORMAT( 1x,///,15x,'ENTER INTEGER NUMBER FOR SELECTION ===> ') 
CALL CUP(21,56) 

READ(*,*) SEL 

IF (SEL .EQ. 1) THEN 


IOPORT = 0 

ELSE IF (SEL .EQ. 2) THEN 
IOPORT = 2 

ELSE IF (SEL .EQ. 3) THEN 
IOPORT = 3 

ELSE IF (SEL .EQ. 4) THEN 
IOPORT = 0 

ELSE IF (SEL .EQ. 5) THEN 
IOPORT = 50 

ELSE 
IOPORT = 99 

ENDIF 


Perecctiln EQ. 4) OR, (SEL .EQ. 5)) THEN 
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CAI CLR 

WRITE(*, 60) 
WRITE(* ,61) 
WRITE(* ,62) 
WRITE(*,61) 
WRITE(*, 63) 
WRITE(* ,61) 
WRITE(*,93) 

60 FORMAT(1x,///,10x,' rer r rrr rer nr rrr rrr ne rte nn nn nn ee nee e--- 


ot lGx ee) 
61 FORMAT( 10x, '-------------2222- 2-222 ee enn n nnn e eee e eee ree econ ene n eee 


) 
62  FORMAT(10x,' ',2x,'SELECT NO. ',3x, ‘,15x, BAUD RATE 9 22=0eeuee 
63 mains Se! 1 b 2k, Vex sec ulOs 
? x 
vo, 10x,” 
* 24x," 
© / dOxae 
2G 
*,/ 10x,’ rep ons 4 (C2 See O00), 
52 ee ne) 
CALL CUP( 18,56) 
READ(*,*) SEL 
IF (SEL .EQ. 1) THEN 
IOPORT = IOPORT + 300 
ELSE IF (SEL .EQ. 2) THEN 
IOPORT = IOPORT + 1200 
ELSE IF (SEL .EQ. 3) THEN 
IOPORT = IOPORT + 4800 
ELSE 
IOPORT = IOPORT + 9600 
ENDIF 


Bore 2 ‘Sox ee Oe 


t 


Soha 3 2x aloo GOO) 


be “ Ye 


CATE CER 
WRITE(*, 50) 
WRITE(**,51) 
WRITE(*,52) 
WRITE(**,51) 
WRITE(* ,53) 
WRITE(**,51) 
WRITE(* ,93) 
50 _FORMAT( 1x, [1] VOX," 92-9908 2 ene 


Oo ') 
Sal FORMAT( 10x, ' De 
52 FORMAT(10x,' ',2x,'SELECT NO. °,3x,  “,17, PARIIY 623e® 
53 FORMAT( 10x, ’ "3X, 1 "2x, |, és NO Pama 

ox 

vo, 10x, ' | ars > ore 2 ',2x,' ',14x, EVEN Parte. 

te 

*,/, 10x, es 5 ony 3 '2xy' * 514%, ORD erases 

J 2 ee) 

CALL CUP( 17,56) 


66 


READ(*,*) SEL 

END IF 

IF (SEL .EQ. 3) THEN 
IOPORT = IOPORT + 1 

ELSE IF (SEL .EQ. 2) THEN 
IOPORT = IOPORT + 2 

EENDIF 

RETURN 

END 


C SUBROUTINE SIZE -- ALLOWS ADJUST OF OUTPUT PLOT * 


15 


16 


ye 


SUBROUTINE SIZE (FACT) 
REAL FACT 
ADJUST SIZE OF OUTPUT PLOT 

CALL CLR 

Wrote, 15) 

FORMAT( 1x, /T/1 10) ne 
Kencene ba / 10x," wie This is the routine to Adjust the plot size deve 
A A pee 
*///,10x, ‘Output plots are currently sized to fill the console scre 
men. ,//, 

* 10x, "Would you like to adjust the plot size?’ ,/,10x, 
vs "Enter Ny" or ‘n' ==> ') 

CALL CUP( 14,31) 

READ(*,'(A)') CHG 

maccGHG .EQ. ‘n’) .OR. (CHG .EQ. ‘N')) THEN 

GQ TO 17 

END IF 

WRITE(*, 16) 

FORMAT(///,10x,'A default plot factor size or 049° 1S jised) to £111 
phe @ensole screen. ,/,10x, 

‘This plot size will be halved by entering a value of Oec5../ , 

*10x, 'It will be doubled with a faCtormemEnya on wl o a 74 10x, 
¥ Venter f4etouisize ——> | ) 

ebb GUP( 22,32) 

READ(*,**) FACT 

RETURN 

END 
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APPENDIX C. D.S.L SIMULATION PROGRAM AND DATA FILE. 


This appendix contains three listing files and its resulting data files for the D.S.L 
simulation programs used in examples 4-1, 4-2 and 4-3. These programs are used to 
produce the data files which generate the plots of the performance index curves on the 


parameter plane. 


A. SIMULATION PROGRAM FOR EXAMPLE 4-1. 


TITLE PERFORMANCE INDEX CALCULATION 

ARRAY Z(1),P(4),C(2),D(2) 

TABLE Z(1)=100. ,P(1-4)=1.0,11.0,10.0,0.0,D(1)=1. 0 

PARAM A=0. 0,B=0. 0, Y0=0. 0 

INITIAL 
6(1)=A 
C(2)=B 
D(2)=B 
BETA=B 
ALPHA=A 

DYNAMIC 
R=STEP(0. ) 

DERIVATIVE 
Y1=TRNFR(0,3,Y0,Z,P,E) 
M=TRNFR(1,1,Y0,C,D,Y1) 
E=R-M 
ABE=ABS(E) 
PI1=INTGRL( YO ,E*E) 
PI2=INTGRL( YO, ABE*TIME) 
PI3=INTGRL( YO, ABE) 

CONTROL FINTIM=15 ,DELT=0. 01 

PRINT 15. , ALPHA, BETA, PI1,PI2,PI3 

END 

PARAM A=4. ,B=10. 

END 

PARAM A=5. ,B=7.5 

END 

PARAM A=5. ,B=10. 

END 

PARAM A=5. ,B=15 

END 

PARAM A=6. ,B=10 

END 

PARAM A=7. ,B=30 

END 

PARAM A=7. ,B=40 

END 

PARAM A=7.5,B=15 

END 

PARAM A=7.5,B=20. 
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END 

PARAM A=7.5,B=25. 
END 

PARAM A=7.5,B=35 
END 

PARAM A=7.5,B=45. 
END 

PARAM A=7.5,B=50. 
END 

PARAM A=10. ,B=10. 
END 

PARAM A=10. ,B=15 
END 

PARAM A=10. ,B=20. 
END 

PARAM A=10. ,B=25. 
END 

PARAM A=10. ,B=35 
END 

PARAM A=10. ,B=40. 
END 

PARAM A=12.5,B=10. 
END 

PARAM A=12.5,B=15 
END 

PARAM A=12.5,B=20. 
END 

PARAM A=12.5,B=25. 
END 

PARAM A=12.5,B=35 
END 

PARAM A=12.5,B=45. 
END 

PARAM A=12.5,B=50. 
END 

PARAM A=15. ,B=10. 
END 

PARAM A=15. ,B=15 
END 

PARAM A=15. ,B=20. 
END 

PARAM A=15. ,B=25. 
END 

PARAM A=15. ,B=30. 
END 

PARAM A=15. ,B=35 
END 

PARAM A=15. ,B=40. 
END 

PARAM A=15. ,B=45 
END 

PARAM A=15. ,B=50. 
END 

PARAM A=17.5,B=15 
END 

PARAM A=17.5,B=20. 


END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 


A=17. 
A=17. 
A=17. 
A=20. 
A=20. 
A=20. 
A=20. 
A=20. 
A=20. 
A=20. 
A=20. 
A=25. 
A=25. 
A=25. 
A=25. 
A=25. 
A=25. 
A=30. 
A=30. 
A=30. 
A=30. 
A=30. 
A=50. 
A=50. 
A=50. 
A=50. 
A=50. 


A=50. 


5, B=25. 
5 ,B=35 


5 ,B=40. 


,»B=15 


, B=20. 


,B=2Z5. 


, B=35 


, B=40. 
, B=30. 


, B=60. 


, B=90 


,»B=15 


, B=20. 


,o—25. 


, B=35 


, B=40. 


, B=25. 


, B=30 


, B=35. 
,B=40. 


,B=50. 
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END 
PARAM 
END 
PARAM 
END 
PARAM 


PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
SLOP 


A=50. 
A=60. 
A=60. 
A=60. 
A=80. 
A=80. 


A=80. 


, B=70. 


A=80. , B=90. 


A=100. 
A=100. 
A=100. 
A=100. 
A=100. 
A=100. 
A=120. 
A=120. 
A=130. 
A=130. 
A=140. 


A=140. 


(ESO: 
, B=40. 
, B=50. 
, B=60. 
, B=80. 
,B=100. 
,B=35. 
BESO. 
,B=35. 
(BESO. 
, B=35. 


,B=50. 


al 


DATA FILE FOR EXAMPLE 4-1. 


el ey 


0.0000 
4.0000 
5.0000 
5.0000 
5.0000 
6.0000 
7.0000 
7.0000 
7.5000 
7.5000 
7.5000 
7.5000 
7.5000 
7.3000 
10.000 
10.000 
10.000 
10.000 
10.000 
10.000 
12.500 
12.500 
S00 
E2500 
12.500 
12.500 
15.000 
15.000 
15.000 
15.000 
15.000 


inst IAS 


0.000 
10.000 
7.500 
10.000 
15.000 
10.000 
30.000 
40.000 
15.000 
20.000 
25.000 
35.000 
45.000 
50.000 
10.000 
135.000 
20.000 
25.000 
35.000 
40.000 
10.000 
15.000 
20.000 
25.000 
45.000 
50.000 
10.000 
15.000 
20.000 
25.000 
30.000 


Les 
1.00000 
0.32856 
0.30910 
0.28481 
0.30855 
0.26177 
O57 25 
0.41387 
0.23340 
0.25469 
0.28477 
Omer, 
0.42382 
0.45935 
0.25001 
0.20401 
0.21075 
OR 22105 
0.27603 
0.30143 
0.28051 
0.19186 
0.18744 
0.19806 
0.27154 
O.2919s 
0.34726 
0.18889 


0.17409 


0.17885 
0.18996 


vave |! 
1.50000 
0.47462 
0.53308 
0.333823 


oS? 
Oe oles 
Orr 3 U7 
0.76717 
1.27460 
Wes) 0) 
0.42101 
Osis 25 
0.14365 
0.19092 
0.38136 
0.50587 
0.66678 
0.16546 
0.11826 
0.12439 
0.38180 
0.47837 
P2710 
0.18998 
0.11261 
0.10424 
0.11504 


Lele 
1.00000 
0.64140 
0.64346 
0.54590 
0.60431 
0.50802 
Cai Weve 
0.92510 
0.42445 
0.47959 
0.56794 
0.76028 
0.95456 
1.05080 
0.54982 
Oy i> 
Or37 509 
0.42358 
0.55835 
UO2725 
0.66477 
Ore 0 
0.33404 
0739079 
0.55436 
0.61059 
0.87889 
0.38491 
0.31716 
0.31487 
0.33643 


15.000 
15.000 
15.000 
15.000 
17.500 
17.500 
17.500 
7-500 
17.500 
20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
20.000 
25.000 
25.000 
25.000 
25.000 
25.000 
25.000 
30.000 
30.000 
30.000 
30.000 
30.000 
50.000 
50.000 
50.000 
50.000 
50.000 
50.000 
50.000 


35.000 
40.000 
45.000 
50.000 
15.000 
20.000 
25.000 
35.000 
40.000 
15.000 
20.000 
25.000 
35.000 
40.000 
30.000 
60.000 
90.000 
15.000 
20.000 
25.000 
35.000 
40.000 
30.000 
20.000 
25.000 
30.000 
35.000 
40.000 
20.000 
25.000 
30.000 
35.000 
40.000 
50.000 
60.000 


0.16644 
0.16621 
0.18434 
0.19669 
U-20112 
0.16250 
0.15766 
0.17018 
0.18022 
0.16203 
0.22763 
0.30424 
23799 
0.16216 
0.14800 
alo 
0.15796 
0.14718 
0.16949 
0.14448 
0.13889 
Ole 70) 
OFT4S9 2 
C1020 
OL NGS UE, 
0.13466 
0.12415 
0.12067 
Onl 200 
le. 


0.14462 
0.19111 
0.24730 
0.31143 
0.23148 
0.11377 
0,09545 
0.10984 
0.13387 
0.29311 
0.12059 
0.09084 
0.09490 
0.10688 
0.08883 
0.23630 
0.56953 
0.51682 
0.15589 
0.09050 
0.07941 
0.08447 
0.07807 
0.21232 
0.10924 
0.07458 
0.06935 
0.07283 
1.10170 
0.26752 
0.15155 
0.10402 
0.07715 
0.05161 
0.05862 


13 


OS) Alles, 
0.41479 
0.45996 
0.50594 
0.41186 
Ore 23) 7) 
0.29500 
0.32730 
Oe) Fl) 
0.4512] 
OS isa7 
O29 
O20 0E2 
O22 
0.28474 
0.44709 
0.65381 
0.57663 
Os See 
OL SME, 
20597 
O27 103 
UO2 76 
7700 
0.28776 
51 
0.24805 
0.25209 
0.81409 
Ose 
OF Oe 
DZ 5509 
Orson 
OZ lod 6 


222710 


60.000 
60.000 
60.000 
80.000 
80.000 
80.000 
80.000 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
120.00 
120.00 
130.00 
130.00 
140.00 
140.00 


30.000 
70.000 
90.000 
30.000 
50.000 
70.000 
90.000 
30.000 
40.000 
50.000 
60.000 
80.000 
100.000 
35.000 
50.000 
35.000 
50.000 
35.000 
50.000 


0.14256 
012563 

0.13784 
0.17846 
0.1084] 

0.10938 
0.11809 
0.27014 
0.12869 
0.10764 
0.10224 
0.10374 
0.11050 
0.20936 
0.11045 
0.25156 
Oa 
0.32030 
Onto 


0.20934 
0.05488 
0.07347 
0.38818 
0.10385 
O:0ono77 
0.05032 
0.88639 
O. 2387 
0.14872 
Ol 17 
0.06032 
0.04206 
0.58018 
O.19035 
O.SO858 
0.22142 
1.26700 
0.24697 


74 


0.338520 
0.21479 
0.24487 
0.45147 
0.22883 
Ome o2 
0.20455 
0.71319 
0.31758 
0.24754 
0.21883 
0.1954] 
OUN2aS0 
0.5493 1 
On26052 
0.66742 
0.28185 
0.85656 
OL29a24 


— 


C. SIMULATION PROGRAM FOR EXAMPLE 4-2. 


TITLE PERFORMANCE INDEX CALCULATION 

ARRAY A(1),B(4) 

TABLE A(1)=50. ,B(1)=1. 0,B(4)=0. 0 

PARAM K1=0. 2,K2=0. 0,Y0=0. 0 

INITIAL 
B( 2)=34+50%K2 
B( 3)=2+50"K1 
BETA=K2 
ALPHA=K1 

DYNAMIC 
R=STEP(0. ) 

DERIVATIVE 
M=TRNFR(0,3,Y0,A,B,E) 
E=R-M 
ABE=ABS(E) 
PI1=INTGRL( YO ,E*E) 
PI2=INTGRL( YO, ABE*TIME) 
PI3=INTGRL( YO, ABE) 

CONTROL FINTIM=15 ,DELT=0. 01 

PRINT 15. ,ALPHA,BETA,PI1,P12,P13 

END 

PARAM K1=0. 2,K2=0. 1 

END 

PARAM K1=0. 2,K2=0. 2 

END 

PARAM K1=0. 2,K2=0. 3 

END 

PARAM K1=0. 2,K2=0. 4 

END 

PARAM K1=0. 2,K2=0.5 

END 

PARAM K1=0. 2,K2=0. 6 

END 

PARAM K1=0. 2,K2=0. 7 

END 

PARAM K1=0. 2,K2=0. 8 

END 

PARAM K1=0. 2,K2=0. 9 

END 

PARAM K1=0. 2,K2=1. 0 

END 

PARAM K1=0. 3,K2=0. 1 

END 

PARAM K1=0. 3,K2=0. 2 

END 

PARAM K1=0. 3,K2=0. 3 

END 

PARAM K1=0. 3,K2=0. 4 

END 

PARAM K1=0. 3,K2=0.5 

END 

PARAM K1=0. 3,K2=0. 6 

END 

PARAM K1=0. 3,K2=0. 7 


1 


END 
PARAM 
END 
PARAM 
END 
PARAM 


PARAM 


K1=0. 
K1=0. 


K1=0. 


K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 
K1=0. 


1 K1=0. 


K1=0. 


K1=0. 


KOS. 
5) =O). 
koe 
oar 
4,,K2=0. 
hn) 
4,K2=0. 
BOC. 
A 0). 
hI). 
4,K2=0. 
=). 
5 ,K2=0. 
5 ,K2=0. 
5 =I). 
5 LO). 
k= 
5 ,K2=0. 
Cy 0. 
LIE, 
ne CC) 
Ps EM 
.6,K2=0. 
.6,K2=0. 
.6,K2=0. 
6 ,K2=0. 
6 ,K2=0. 


.6,K2=0. 
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K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 
K1=0 


K1=0 


.6,K2=0. 
.6,K2=0. 
.6,K2=0. 
non koa 
ey ko=o. 
By kO=0: 
a7 K2=0, 
TAG 
k=O. 
ko 
ek o=0. 
yk 2—0: 
7, K2=0. 
7 ,K2=1. 
8, K2=0. 
Boat. 
2 TO). 
.8,K2=0. 
8, K2=0. 
oko On 
wemk2—OR 
none —or 
.8,K2=0. 
oily 
non2—or 
OK). 
Wo k2=0) 


.9,K2=0. 


nn F&F WW 


ne 


END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 


SLOP 


l=). 6) K=O. 
G0), 8), P—-9)- 
0), ©), P10). 
K1=0. 9 ,K2=0. 
MEO, O° 0) 
k=O eoeo—de 
lik, 0,210) 
GSE), 0) PSO. 
EM, 0) 2S. 
Kil=1)0nk—0; 
Ki= Tonk o=o- 
ily ©) KP. 
Ki loeKe 
le KOS: 
as (0), RSC: 
K1=1. 0,K2=0. 


K1=1. 6,K2=1. 
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Peer riA 


0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.20000 
0.30000 
0.30000 
0.30000 
0.30000 
0.30000 
0.30000 
0.30000 
0.30000 
0.30000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.40000 
0.50000 


ajc jan 


0.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1 .Q0000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.S0000 
0.90000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 
0.10000 


D. DATA FILE FOR EXAMPLE 4-2. 


ES 
1823.40 
0.81565 
0.91714 
1.09550 
1.28850 
1.48580 
1.67890 
1.87450 
2.06100 
P25 130 
2.40750 
0.54210 
0.66416 
0.80281 
0.94561 
1.09000 
1.23450 
1.37940 
132230 
1.66200 
0.47397 
0.57805 
0.68821 
0.80004 
0.91257 
1.02540 
1.13840 
1.25100 
le36320 
1.47550 
0.46277 


PS vose 


1177.4 
4.9245 
6.0067 
8.5450 
11.4640 
14.6030 
17S 0 
20.270 
Balen gal, 
24.4710 
26.7370 
5228 
2.3429 
3.8884 
05 
S10 
PO SO0 
f23 150 
14.3720 
16.1900 
0.5706 
ree | 
ZA So) 
3.3445 
4.7404 
7507 
7.9484 
0S! 
11.1670 
13.0080 
0.3224 


te 


0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.60000 
0.70000 
0.70000 
0.70000 
0.70000 
0. 70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.$0000 
0.80000 
0.80000 
0.80000 
0.80000 
0.80000 


0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.$0000 
0.90000 
1.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 
0.10000 
Q. 20000 
0.30000 
0.40000 
0.50000 
0,60000 


Ois5003 
0.64156 
On so22 
0.82524 
0.91744 
1.00970 
1.10200 
1.19420 
1.28650 
0.47534 
0.55087 
0.62798 
0 US 
0.78336 
OENGilS 
0.93921 
EO1720 
1.09520 
i bias 20) 
0.50008 
0.56605 
Ores 29% 
0.70021 
OPS ccwieyy 
0.83499 
0.90245 
OGY 3 
1.03740 
1.10490 
Ose)? 
OFSI030 
0.64946 
0.70875 
0.76813 
0.82756 


Os 2 
1.3845 
Zales) 
Solo? 
4.2318 
5.4501 
6.7009 
8.0218 
9.4955 
WS 
0.5490 
Ores 
ES o79 
2.2462 
3.0514 
SE 525 
4.9258 
Dot 2 
fA ae 
0.3878 
0.4620 
0.7805 
1.2014 
bag bola) 
2.5256 
3.0156 
3.7809 
4.6278 
Spo 
0.5459 
0.4757 
Oro7 91 
ORS eis 
eo9278 
1.8649 
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0.9005 
le) 
eer o2 
e610 
1.8455 
2.0760 
2.2984 
2.5160 
2.17345 
0.6421 
0.8301 
BO 257 
eZ 
1.4262 
1.6245 
ol 
2.0147 
2 0Gr 
Peis Mis 
0.7400 
0.8138 
0.97306 
1.1408 
1.3100 
1.4801 
1.6500 
1.8192 
1.9882 
Wa) SS ve 
0.8400 
0.8496 
Ws SoM 
1.0967 
1.2406 
Looe 


0.80000 
0.80000 
0.80000 
0.80000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
0.90000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.6000 


0.70000 
0.80000 
0.90000 
1.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 
0.10000 
0.20000 
0.30000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 


0.88701 
0.94648 
1.00600 
1.06540 
0.56816 
0.62062 
OO sv 
O05 5 
0.77964 
0.83276 
0.88590 
O75 205 
W922 1 
1.04540 
0.60743 
0.65498 
0.70284 
0.70284 
0.75080 
0.79880 
0.84683 
0.89487 
0.94292 
OPO O09 7 
1.14690 


2.4071 
5.0174 
5.6918 
4.414] 
On259 
0.6242 
0.6623 
0.8877 
1.1955 
1.5688 
2.0029 
2.4950 
3.0896 
SE0208 
OP) ie 
0.8218 
0.7427 
0.7427 
0.8598 
1.0886 
1.3838 
1.7344 
21352 
2.5824 
2 Oren 
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| ops)oys 
1.6838 
1.8324 
1.9801 
0.9400 
0.9400 
Vie2s 
1.0862 
Ie AO)S) 
ee 
1.4606 
eo 
PAS 
1.8540 
1.0400 
1.0400 
1.0448 
1.0448 
1.1062 
feleee 
05> 
1.4169 
Poole 
1.6476 
1.7147 


E. SIMULATION PROGRAM FOR EXAMPLE 4-3. 


TITLE PERFORMANCE INDEX CALCULATION 
ARRAY A(1),B(3),C(2),D(2) 
TABLE A(1)=0. 1084,B(1-3}=90. 36,1.0,0. 0,DC1)=1. 0 
PARAM Z=0.0,P=0.0,Y0=0.0 
CONST K1=0. 454616, R=16. 296 
INITIAL 
C(1)=Z 
C(2)=P 
D(2)=P 
BETA=P 
ALPHA=Z 
DYNAMIC 
R=STEP(0. ) 
DERIVATIVE 
SI=TRNFR(0,2,Y0,A,B,DEL) 
DEL=TRNFR(1,1,Y0,C,D,X1) 
SIE=SI-R 
X1=K1*"SIE 
SIES=SIE*“SIE 
DELS=DEL*DEL 
=R*‘SIES+DELS 
J=INTGRL(Y0,Q) 
CONTROL FINTIM=15. ,DELT=0. 05 
PRINT 15. ,ALPHA,BETA,J 
END 
PARAM Z=2. ,P=0. 1 
END 
PARAM Z=2. ,P=0. 3 
END 
PARAM Z 
END 
PARAM Z=2. ,P=0. 7 
END 
PARAM Z=2. ,P=0.9 
END 
PARAM Z=2. ,P=1.0 
END 
PARAM Z=4.5,P=0. 1 
END 
PARAM Z=4.5,P=0. 3 
END 
PARAM 2=4.5,P=0.5 
END 
PARAM Z=4.5,P=0.7 
END 
PARAM Z=4.5,P=0.9 


2. ,P=0.5 


PARAM Z=4.5,P=1.0 
PARAM Z=7.5,P=0. 1 
PARAM Z=7.5,P=0. 3 


PARAM Z=7.5,P=0.5 


END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 


PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 


Z=7. 
Z=/. 
Z=7. 
Z2=9. 
Z=9. 


Z=9. 


Z=9. 
Z=9. 
Z=13. 
Z=13. 
Z=13. 
Z=13. 


Z=13. 


Z=19. 
Z=19. 
Z=19. 


Z=23. 


One 
5,P=0. 9 


Sigal 


.5,P=0. 
.5,P=0. 


.5,P=1. 


END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 


STOP 


Z=23. 
Z=23. 
Z=23. 
Z=23. 
Z=23. 
Z=26. 
Z=26. 
Z=26. 
Z=26. 
Z=26. 
Z2=26. 


Z=29. 


Z=29. 
Z=29. 
Z=35. 
Z=35. 
Z=35. 
Z=35. 
Z=35. 


Z=35. 


Bee). 3 
5 ,P=0. 5 
5) Se, 7 
5 ,P=0.9 
5 Pi 
ol 
,P=0. 3 
pP=0e5 


,P=0.7 
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ALPHA 


0.0000 
2.0000 
2.0000 
2.0000 
2.0000 
2.0000 
2.0000 
4.5000 
4.5000 
4.5000 
4.5000 
4.5000 
4.5000 
7.5000 
7.5000 
7.5000 
7.5000 
7.5000 
9.0000 
9.0000 
9.0000 
9.0000 
9.0000 
9.0000 
13.0000 
13.0000 
13.0000 
13.0000 
13.0000 
13.0000 
16.5000 


BELA 


0.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 


F. DATA FILE FOR EXAMPLE 4-3. 


J 

15.000 
14.937 
15.050 
15.042 
e020 
15.002 
14.994 
14.899 
14.982 
14.982 
14.970 
14.959 
ole a) 
14.862 
14.912 
14.918 

14.912 
14.91] 

14.846 
14.881 

14.889 
14.890 
14.890 
14.890 
14.812 
14.809 
14.820 
14.829 
14.836 
14.839 
14.79] 
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16.5000 
16.5000 
16.5000 
16.5000 
16.5000 
19.0000 
19.0000 
19.000 
19.000 
19.000 
19.000 


nr oO tv td Ft 
ow tos Us td Gs or 
Cry Um 
Co aa, er SS 
a, CS pe ES Fae EP a OE 


ho 
(sy Un 
OS 
=) 


26.000 
26.000 
26.000 
26.000 
26.000 
26.000 
29.500 
29,900 
29.500 
29.500 
272200 
29.500 
35.000 
35.000 
35.000 
35.000 
35.000 
35.000 


0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
0.10000 
0.30000 
0.50000 
0.70000 
0.90000 
1.00000 
1.00000 
0.30000 
0.50000 
0.70000 
0.90000 
1 .QOOOO 


14.759 
14.769 
i782 
14.794 
14.799 
14.779 
14.730 
14.737 
14.752 
14.766 
Bone 
14.765 
14.686 
14.689 
14.704 
14.722 
14.730 
14.760 
14.667 
14.666 
14.681 
14.699 
14.708 
14.755 
14.646 
14.638 
14.652 
14.671 
14.680 
14.754 
14.621 
14.603 


14.614 


14.632 
14.642 
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APPENDIX D. USER’S MANUAL 


When we wish to plot the constant cost curves, the following described procedures 


should be performed as indicated. In this section, let us consider the specific problem 
in example 4-1. 


fmeeeolliinnG [itE PARAMETER PLANE CURVES. 


The first step is to plot the parameter plane curves which are constant zeta curves 
and constant omega curves as shown in Figure 12 in Chapter 4. This curve permits us 


to choose the range of two parameters that are to be used in the analysis and design. 
PeeteNERATING A DATA FILE. 


The second 1s to obtain a data file containing a , 8 and performance index values. 
Simulation in D.S.L is used to evaluate the cost at a large number of points on the pa- 


rameter plane. The sinulation results look as follows. 


1 DOU/VS ....4 RELEASE 1 ..... MODIFICATION LEVEL 1 ..... 1985 
0 
0 >>>> DSL SIMULATION INPUT DATA <<<< 


TITLE PERFORMANCE INDEX CALCULATIOM 

TABLE 201)=100.,P01-4)=1.0,11.0,10.0,0.0,D€1)21.0 
PARAM A=0.0,B=0.0,Y0=0.0 

CONTROL FINTIN=A15, DELT=0.01 

PRIWE 15.,ALPUA, DETA,1ES, 1AET,1AE 


Ettp 
OxKN KCLENGH = 2189, KPOINME = 227, AVAILABLE SPACE LEFT IN COMMON/CURVALZ = 
OxXxX RKS IMTEGRATION METHOD USEN xxx 
0 REX DSL OUIPUT LISTING, GROUP 1 
+ » STARTING RUN 1 KKK 


PERFORMANCE IMDEX CALCULATION 
ALPHA BETA 


TINE IES IAET IAE 
0.00000E#00 O.0000DE!00 O.ODD0D0D0EF00 O.C0000E+00 O.00000Et00 0,00000E+00 
, 15.000 0.00000Et00 0.00000E#00 15.000 112.50 15.000 
0 >>>> DSL SIMULATION LMPUT DATA <<<< 
PARAW A=4.,D210. 
Ettb 
0 *x¥M DSL OUIPUT LISTING, GROUP 1 
- » STARTING RUIN 2 KR 
PERFORMANCE INDEX CALCULATION 
TIVE ALPHA BEIA IES LAET LAE 
0.00000E+00 4.0000 30.000 0.00000E#00 O.00000Et00 o0.00000E+00 
. 15.000 4.0000 10.000 0.32856 0.974962 0.649140 


We have to edit the data file as shown in Appendix C to remove unnecessary parts. 
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3, LOADING AW Aa 


The third step is to download this file from the mainframe computer to a floppy 
diskette. This data now on the floppy diskette can be copied on the IBM microcomputer 
under a designated filename. In this example, the filename is called pil.dat. Now, we 
are prepared to plot the families of performance index curves in the same range of pa- 


rameter values on the two parameter plane. 
4. PLOTTING THE COST Guia es 


To run the performance mdex program, type in the executable file name, 
COSTFUN. The IBM microcomputers in the NPS Controls Lab are set up under the 
file management system IDIR. To execute the program under the this system, enter the 
appropriate subdirectory by using the up/down arrow keys to position the cursor adja- 
cent to the subdirectory name, then depress the ENTER kev. Use this same procedure 
to select COSTFUN.EXE and push ENTER to run the program. The preparation step 


for plotting the desired curves will appear. 


Enter proper values after each prompt as follows: 


1.llow many constant PI curves # ===>6 

2.fnter PI values ,lwgt,ldig # 12===20.27,1,2 
2.Enter PI values ,lwgt,ldig # 2===>0.2,1,1 
2.Enter PI values ,lwgt,ldig # 3===50.15,1,2 
2.Enter PI values ,lwgt,ldig # 4===590.12,1,2 
2.Enter PI values ,lwgt,ldig # 5===>0.11,1,2 
2.Enter PI values ,lwgt,ldig # 6===>0.1,1,1 
3.Enter number of data file points ===>86 
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Next, adjust the plot size to the same parameter plane curves. 


Output plots are currently sized to fill the console screen. 


Would you like to adjust the plot size? 


"1 | " ? 


Enter "yor on ==>) y 


A default plot factor size of 0.9 is used to fill the console screen. 
This plot size will be halved by entering a value of 0.45. 
It will be doubled with a factor entry of 1.8. 


Enter’factor size ==> 0.7 


The last step is to enter a data [tle name. 


File name should not exceed 8 characters 
File extension should not exceed 3 characters 


What is the file name (fn) and extension (ext)? 


Enter in form fn.ext (e.g. PARAPLAN.INP) ===> pil.dat 
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After displaying the graph on the monitor, the program prompts you for a hardcopy 
of the graph. 


WOULD YOU LIKE HARDCOPY OF THESE PLOTS? (Y or N) 


-— ee ee ew om ow ow ae ame oe ee om ee ee om oe eo oo om ee em ee em em 6 ew em om ee om em ew ew of ee @ =8 oF ss of oe of ee of ee of eo ef of ew =8 of ew ef == of oe ow ow «o8 «8 == oe = 


; 1 Epson FX-80, All 4 
: 2 ' Epson FX-100, All : 
: 3 Epson MX-100, All ; 
4 : Epson RX-80, All ; 
. 5 Epson MX-80 & IBM Printer 
: 6 HP 7470A Graphics Plotter : 
: il } HP 7475A Graphics Plotter : 
4 8 : HP 758xB Series Plotters 
: 9 : HP 2686A Laser Jet : 
HEFEHEFEFEFE EEE HEE EEE EEL EEE EEE FEE HEF EF EEF FEF HEHE EEE EEE HH Hee t tet 
: 10 : Graphics Monitor (default) ; 
: 11 ; HARDWARE Interface Menu ' 
' 12 : Input PLOT88 Values for IOPORT and MODEL jj; 
! 99 : RETURN to LOAD Data : 


Enter integer number for selection ===> 5 


Reenter the data file name for the hard copy. Figure 13 is a performance index 
curve for constraint f£? dt. In like manner, the other two nonlinear criterion of per- 
formance index in terms of f1E| ¢dt and flE| dt can also be plotted by running the 
same program. Figure !4 and Figure 15 are shown for desired constraint { |E| «dt and 
{lL dt, respectively. 
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